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Abstract 

An explicit expression in terms of canonical variables is obtained for the Hamiltonian functional 
determining the fully nonlinear dynamics of two-dimensional potential flows of an ideal fluid with 
a free surface over an arbitrary nonuniform depth. The canonically conjugate variables are derived 
from the previously developed non-canonical conformal description of water waves over a strongly 
undulating bottom [V. P. Ruban, Phys. Rev. E 70, 066302 (2004)]. Also an alternative approach 
to the problem is discussed, which gives weakly nonlinear Hamiltonian models of different orders. 
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I. INTRODUCTION 



For many phenomena related to the science of ocean waves, it is necessary to study fully 
nonlinear regimes in the dynamics of free water surface. For instance, rogue waves are 
highly nonlinear structures attracting much interest over past years Breaking waves 
on a beach give another example of nonlinear water wave dynamics. For some problems 
the simplified mathematical formulation can be used, when the fluid is treated as inviscid 
and incompressible, and only potential flows with a constant pressure above the free surface 
are considered. Such class of ideal flows is known to possess the Hamiltonian structure j^J, 
where the pair of canonically conjugate functions is formed, for example, by the vertical 
displacement 77 (r, t) of the free surface from the mean horizontal level and by the boundary 
distribution ip(r,t) of the velocity potential, with r being the position in the horizontal 
plane and t being the time. Unfortunately, the variables r\ and ip are not suitable for highly 
nonlinear waves, since there is no explicit expression for the Hamiltonian functional in terms 
of 77 and ip, valid for arbitrary wave steepness. For two-dimensional (2D) flows (in xy-plane) 
this difficulty can be avoided, and exact and compact equations of motion are possible with 
the help of so called conformal variables originally introduced by L.V. Ovsyannikov in early 
1970's 3|,y|. The conformal description was combined with a variational formalism in 1990's 
by V. E. Zakharov and co-workers, for the infinite depth {5] and for the case of flat horizontal 
bed Subsequent development of this method was done in Refs.[?-10]. The general idea 
of the approach is to use a time-dependent conformal mapping x + iy = z(w, t) which 
"straightens" the flow domain, with an analytic function z(w,t) depending on a complex 
variable w = u + iv . The free surface and the bottom correspond to some fixed values v s and 
Vb of the curvilinear coordinate v (in the case of infinite depth v s — and Vb = —00, while in 
the case of finite depth without loss of generality one can take v s = 1 and Vb = 0), and thus 
the conformally invariant 2D Laplace equation for the velocity potential is easily solved. 
Basic unknown dynamical quantities, such as the boundary value ip of the velocity potential 
and the vertical coordinate Y along the free surface, in this description are functions of u 
and t. Explicit expressions for time derivatives ip t (u,t) and Y t (u,t) can be derived from 
variational Euler-Lagrange equations with the help of a linear integral relation between 
Y(u,t) and X(u,t), taking place due to the analyticity [5|, |6|. 

Non-trivial bottom topography is known to affect wave motion. In order to include into 
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the fully nonlinear model that influence, the conformal description was generalized by the 



present author to the case of static nonuniform bed ll], for time- dependent bottom bound- 
ary [l^, and for planar flows with constant vorticity over nonuniform bed 13]. Even for 
weakly nonlinear regimes, conformal coordinates in the problems with bottom topography 
are more preferable than the Cartesian coordinates, though in that case a somewhat different 



approach is employed (see, for example, Refs 



15], and Sec. Ill below). A step to- 



wards using conformal variables for three-dimensional (3D) highly nonlinear potential flows 



with a free surface was done in the paper 



161 ]. where for long-crested waves an asymptotic 



expansion of the Hamiltonian functional was suggested in the powers of a small parameter 
corresponding to the squared ratio of a typical wave length A to a long typical distance l q 
along wave crests stretching in the transversal horizontal direction q. 

Equations of motion in the conformal description contain integral operators which are di- 
agonal in the Fourier representation. Therefore in the numerical simulations a Fast Fourier 



Transform (FFT) can be used (see, for example, the web site [17[ about a modern and 
efficient FFT library). Important numerical results concerning planar rogue waves on in- 
finitely deep water have been obtained by this method in Refs. 18l-l2l|. Nonlinear waves over 
inhomogeneous static and time-dependent beds were simulated in the conformal variables 



in Refs. 



11 



121 ] . The most recent application of this method can be found in Refs. 22|, |23 | 



where the phenomenon of so called water-waye gap solitons over spatially periodic beds was 
studied. The weakly 3D conformal theory 16] has been recently successfully applied to 



24 



27]. 



simulate long-crested freak waves at infinitely deep water 

It should be noted that in most cited papers non-canonical conformal variables were used, 
as they are quite sufficient for numerical simulations. The exceptions are the analytical works 
js, 9], where a pair of canonical conformal variables has been employed, however in the case 
of infinite depth only. For uneven beds, the canonical pair was not derived until now. In the 
present work, this gap in the theory will be filled and canonical conformal variables together 
with an explicit Hamiltonian will be suggested for 2D ideal potential flows over an arbitrary 
nonuniform bottom profile (in Sec. II). But we will see the obtained exact Hamiltonian con- 
tains a strongly non-local linear operator which at long scales is similar to the integrating 
operator d~ Y . The presence of such operator makes the exact canonical model somewhat 
difficult for analysis. Therefore an alternative approach how to introduce canonical variables 
for waves over strongly undulating bottom will be also considered (in Sec. Ill), which is free 
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of the above mentioned technical difficulty. Unfortunately, the alternative approach does not 
allow us to obtain an exact Hamiltonian, but approximations of different orders for weakly 
nonlinear waves are possible. In this paper, the calculations will be made up to the fifth 
order. 



II. EXACT HAMILTONIAN FOR WAVES OVER NONUNIFORM DEPTH 

It is appropriate to mention here that if the given bottom boundary is non-flat, then 
the conformal mapping z(w, t) can be represented as the composition z(w, t) = Z(((w, t),t), 
where Z((, t) is a known analytic function [Z(Q does not depend on time if the bed is static, 
which case is assumed in this work]. The intermediate unknown analytic function ((w,t) 
takes purely real values 7(w,£) at the real axis: ((u + Oi,t) = 7(u,t), and / y u (u,t) > 0. 
As a result, the bed profile is given in the parametric form + iY^ = ^(7), where 
— 00 < 7 < +00, while the shape of free surface is determined parametrically by the formula 

X [s) + iY [s] = t))? ^ f) = + hf y (1) 

Since ((u + i,t) is the analytic continuation of a real function 7(w,t) from the real axis, 
a linear integral relation takes place between the real and the imaginary parts of function 
£(u,t): 

C{u, t) = e- k j{u, t) = (1 + iR)p(u, t), (2) 

where k = —id u is the differential operator, e~ k is the operator making the analytic con- 
tinuation, p(u,t) = [cosh k]j(u, t) is another unknown real function, and R = i[tanhk] is a 
linear antisymmetric operator which is diagonal in the Fourier representation: 

Rp(u,t) = j i[t a nh k] Pk (t)e iku — , (3) 

with Pk{t) = j p(u, t)e~ lku du being the Fourier transform of the function p(u, t). The inverse 
operator is T = R^ 1 = — i[cothk}. For long waves, when \k\ 1, operator R is similar 
to d u , while T is similar to integrating operator d~ x . Explicit formulas for the kernels of 
operators R and T in w-representation are given, for instance, in Ref. 111 ]. 

A simple way to obtain the required canonical pair is the following. We first recall the 
expression for the Lagrangian functional describing water waves over a static nonuniform 
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bed in terms of the conformal variables (see Ref . rjj ) , 



C = / ^\Z\i)\\ Pu R Pt - Pt R Pu ) du - H. 



(4) 



Here Z'(C,) = dZ/d^, and the Hamiltonian "H is the total energy of the system — the kinetic 
energy plus the potential energy in the vertical gravitational field g (surface-tension effects 
can be also included, but we do not consider them in this paper), 

1 



H = - J1>kj> du+ 9 -\ [Im Z(Of Re du, 



(5) 



where K = k tanh k is a Hermitian operator. 

Let us now choose function p(u, t) as the generalized canonical coordinate. Using the 
antisymmetric property of operator R, we rewrite Eq.(j4]) as follows, 



C 



-R 



^\Z'(0\ 2 Pu] - ^\Z\i)\ 2 R Pu ] Pt du - H. 



(6) 



It is clear that the Lagrangian will take the canonical form C = J fipt du — % if we define 
the canonical momentum p by the relation written below 



p = -R MZ\i)\ 2 p u -^\Z\i)\ 2 Rp 



(7) 



Using this equation, it is now necessary to express ip through the canonical variables p and 
p and then substitute the result into Eq.([5]). Generally speaking, this is not a trivial task in 
view of the presence of integral operator R. Fortunately, the coefficients of the equation are 
very special, and therefore the solution can be found in exact form. To solve the integral 
equation ((7j) with respect to ip, we rewrite it in the form 



Re{[l + iR](-iif>\Z'{Z)\ 2 Z u + p)} = 0. 



(8) 



The above equation allows us to conclude that 



i^\z'{i)\% + p = -i{i-ik)f, 



(9) 



where / is some real function. Since £ u = (1 — iR)p u (the overline means the complex 
conjugate quantity), we can multiply Eq.Q by and use a general formula 



[{l-iR)h][{l-iR)f 2 ] = {l-iR)h 



(10) 
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where fx, f 2 , and / 3 are real functions, and / 3 = fif 2 — (Rfi)(Rf 2 ). As the result, we have 



4\Z'(Z)£ v \ 2 + iiii v = (l-iR)f, 



(11) 



where / is another real function. Taking the imaginary part of Eq. ffTTT) . we obtain pp u 
—Rf, and thus / = —T{pp u ). Then the real part gives us the required formula 

-f(pp u ) - pRp u 



(12) 



\d n Z{[l + iR]p)\* 

Now we are able to write explicit and exact expression for the Hamiltonian functional in 
terms of the canonical variables p and p, 



H = 2 
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T {p Pu) + P Rpu 



K 



\d u Z([l + iR}p)\*\ [\d u Z([l + iR]p)\ 



T(pp u ) + pRp u 



du 



+ |/[Im Z([l + iR]p)YRe[d u Z([l + iR]p)]du, 



(13) 



which is the central result of this paper. The obtained formula is rather cumbersome, and 
it is unclear yet if some canonical transformation will be able to reduce this Hamiltonian to 
a simpler form. Further serious work on the problem how to simplify the Hamiltonian ( TT3]) 
is needed. 

The corresponding canonical equations of motion are p t = 5%/8p and —p t = 8%/ 5 p, 
where the variational derivatives should be calculated from Eq.( TT3"j) according to the well- 
established general rules. Let us introduce a short-hand notation 

T(ppu) + pRpu 



N=\d u Z(OrK 



(14) 



Then the variational derivatives are 



— = (Rp u )N - Pu TN, 
op 



(15) 



51-1 

— = d u [R(pN) + pTN] 
dp 



- 2Re 
+ 2Re 
+ g Im 



l-iR){[T{pp u )+pRp u }NZ"{i)/Z'{i)} 
1 - iR)d u {[f{pp u ) +pRp u }N/i u } 

(l-iR){lm(Z(mZ'(0\%} 



(16) 
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III. ALTERNATIVE APPROACH 



An essential difficulty of the derived fully nonlinear canonical model is the presence of 
strongly non-local operator T. Therefore it makes sense to consider alternative approaches 
to the problem of canonical variables for waves over uneven bed. Below we shall focus on the 



approximate method suggested in Ref.[ll|, where curvilinear conformal coordinates u and 
v are static. In this approach, the conformal mapping "straightens" the bottom boundary 
but not the free surface, in contrast to the exact description used in the previous section. 
In Ref. ll|, the corresponding Hamiltonian functional was derived up to the third order in 
canonical variables. Here we develop this method more systematically and calculate the 
fourth- and fifth-order approximations. 

We now introduce (static) curvilinear coordinates u(x, y) and v (x, y), where real function 
v(x,y) obeys the Laplace equation v xx + v yy = with the boundary conditions v = —1 
at the given arbitrary nonuniform bottom and v — at y — 0. Function u(x,y) is taken 
harmonically conjugate for v(x,y). The inverse conformal mapping x + iy = z(u + iv) 
satisfies the condition Im z(u + 0i) = 0, so a real function x(u) = z(u + 0i) parametrizes the 
unperturbed free surface y = 0. The bed profile is given in the parametric form X^+iY^ = 
z(u—i), while the shape of the free boundary is determined through an unknown real function 
V(u,t) by the complex equality 

X [s] +iY [s] = z{u + iV{u,t)). (17) 

The free-boundary value of the velocity potential is ipiu^t) = ip(u,V(u,t),t), where the 
potential tp(u,v,t) obeys 2D Laplace equation ip uu + (p m = in the flow domain —l<v< 
V(u, t), with the bottom boundary condition (p v (u, —1, t) = 0. Consequently, a general form 
of the Lagrangian functional in this description is the following, 

C = J i)\z'{u + iV)\ 2 V t du-H{V^}. (18) 

The Hamiltonian functional 7i{V,ip} = JC{V,ip} + V{V}, where )C{V,ip} is the kinetic 
energy and V{V} is the potential energy. Functional V{V} is relatively simple, 



V = -/[Im z{u + iV)fRe[d u z{u + iV)]du 



V 2 x' 3 + V q 



(x ,2 x"Y 5 



-x' 2 x"' 



du 



+ 0(V 6 ), (19) 



with x' = dx(u)/du, x" = d 2 x(u)/du 2 and so on. As to the kinetic energy JC{V, ip}, it cannot 
be represented in an exact form, but its expansion K, = Kffl + /C^ 3 ' + /C^ 4 ' + ■ ■ • in powers of 
(supposedly small) functions V and if) is possible. We note that the equations determining 
the 2D velocity potential <p(u,v,t) are formally identical to the equations for the velocity 
potential tp(x,y,t) in Cartesian coordinates over a straight horizontal bottom at y — — 1. 
The kinetic energy is determined by the same integral in both cases, 

W 1>) = \j du f^\& + d)dv. (20) 

In view of this equivalence, the expansion of K.{V,ip} is calculated in the same standard 
manner as for weakly nonlinear waves over flat bottom in Cartesian coordinates (compare 
to expansion of fC{rj, in Appendix B of Ref. [231]): 



2 / iP K Hu, 

\ J v[{i, u f - (k^f]du, 
i 

2 



ipkvkvki) + v 2 {kip)ii) v 



du. 



(21) 
(22) 
(23) 



V 3 * * 

—(KtP)k^ 

o 
V 



%l)KVKVKVK%l) 



Yl(k^)dKvkij) 



du. 



(24) 



As a consequence, the bed nonuniformity does not affect ]C{V,ip} in any order. This is an 
essential technical advantage of the conformal coordinates in this problem. 

The form of Lagrangian (|T8|) allows us to treat ip{u) as the canonical momentum if the 
corresponding canonical coordinate xi u ) is defined in the following way, 



xw 



\z'(u + iv)\ dv. 



(25) 



Since we consider small values V, the above equation gives us an expansion of x m °dd 
powers of V, 

x = X ' 2 V + \(x") 2 - x'x"']V 3 /3 + 0(V 5 ). (26) 
From here we can express the non-canonical variable V through the canonical variable %: 



V 



x 



X_ 

a 



3C 30 



[x") 2 \ X 2 



x 



16 



+ o( x 4 



(27) 



S 



Now we substitute the above expression into r H{V,ip} and collect there terms of the same 

order. After simplification, we obtain the Hamiltonian functional %{x, ^} = +"H^ H 

up to the fifth order: 



n [2] 



+ 



9_ 
12 



x' 



du, 



X 



^k^-k^-k^ + ^-ik^ 



X 

2^ 



X 



X 



du 



X 



n\2 



X 



a 



X A du. 



(28) 
(29) 



(30) 



[X 



t/\2 



X 



78 



X 3 du 



+ 



12 



,3 



X 

X 



/6 



X 



X 



x'- 



X 



du. 



(31) 



x'* x'- 

We do not write here the corresponding canonical equations of motion — though it is a 
simple exercise to calculate the variational derivatives 51-1/ 5\ and SH/dip, but the formulas 
are quite long. This weakly nonlinear model can be applied to a wide variety of problems, for 
instance, to study the dynamics of solitons propagating along a channel of variable depth, 
or the scattering of a wave train by a strong bed nonuniformity. The obvious advantages of 
the model are that it is fully dispersive and there is no principal limitations on the bottom 
profiles. However, it is not clear at the moment what is the largest value of wave amplitude 
which is still well described by this approximate Hamiltonian. The answer perhaps depends 
on many details. This question requires further investigations. 

Weakly dispersive (shallow-water) regime can take place if the depth is slowly varying, 
\x"\/x' <^ 1. In that regime one can reduce the Hamiltonian to a purely local form by using 
the long- wave expansion for the integral operator K = —d^ — d*/3 + ■ ■ ■ and then apply 
some analytical methods. Otherwise, the fully dispersive equations of motion can be easily 
simulated on computer with the help of FFT routines, since all linear operators are diagonal 
in Fourier representation. 



9 



To conclude, in this paper two different methods have been presented how to solve the 
theoretical problem of canonical Hamiltonian description of nonlinear planar water waves 
over nonuniform depth. The first method is fully nonlinear, but the corresponding Hamil- 
tonian functional contains a strongly nonlocal operator. The second way is the further 

I I 

development of the weakly nonlinear approach suggested in Ref . . 
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